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ABSTRACT 

We use a simple physical model to study the nonlinear behaviour of the r-mode insta- 
bility. We assume that r-modes (Rossby waves) are excited in a thin spherical shell of 
rotating incompressible fluid. For this case, exact Rossby wave solutions of arbitrary 
amplitude are known. We find that: 

(a) These nonlinear Rossby waves carry ZERO physical angular momentum and 
positive physical energy, which is contrary to the folklore belief that the r-mode angular 
momentum and energy are negative. We think that the origin of the confusion lies in 
the difference between physical and canonical quantities. 

(b) Within our model, we confirm the differential drift reported by Rezzolla, Lamb 
and Shapiro (1999). 

Radiation reaction is introduced into the model by assuming that the fluid is elec- 
trically charged; r-modes are coupled to electromagnetic radiation through current 
(magnetic) multipole moments. We study the coupled equations of charged fluid and 
Maxwell field dynamics and find that: 

(c) To linear order in the mode amplitude, r-modes are subject to the CFS insta- 
bility, as expected. 

(d) Radiation reaction decreases the angular velocity of the shell and causes dif- 
ferential rotation (which is distinct from but similar in magnitude to the differential 
drift reported by Rezzolla et al.) prior to saturation of the r-mode growth. This is 
contrary to the phenomenological treatments to date, which assumed that, prior to 
the saturation of the r-mode amplitude, the loss of stellar angular momentum is ac- 
counted for by the r-mode growth. This establishes, for the first time, that radiation 
reaction leads not only to overall loss of angular momentum, but also to differential 
rotation. 

(e) We show that for I — 2 r-mode electromagnetic radiation reaction is equivalent 
to gravitational radiation reaction in the lowest post-Newtonian order. Based on our 
electromagnetic calculations, we conclude that inertial frame dragging, both from the 
background rotation and from the r-mode itself, will modify the r-mode frequency by 
a factor ~ -Rschwaizschiid/-Rstar7 in qualitative agreement with Kojima (1998). 



1 INTRODUCTION 

Andersson (1998) has shown, and Friedman and Morsink 
(1998) have confirmed analytically, that r-modes of rotating 
stars can grow because of gravitational radiation reaction. 
Lindblom et al. (1998) have shown that this instability can 
be important in rapidly rotating hot neutron stars, where 
the r-mode amplitude might become large enough to af- 
fect the spin frequency [see also Andersson et al. (1998)]. 
The details of nonlinear evolution, which allows for large 
r-mode amplitudes, are essential for astrophysical applica- 
tions [Owen et al. (1998), Spruit (1999), Levin (1999)]; yet 
only phenomenological treatments of such nonlinear evolu- 
tion exist so far. This paper is an attempt to learn about 
the nonhnear behaviour of rotating fluid in which r-modes 
are driven by radiation reaction. 

To address issues of principle, we choose to study a 



very simple system. Our "star" is a thin rotating shell of in- 
compressible inviscid fluid which is sandwiched between two 
hard spheres. These spheres exert no friction on the fluid; 
their role is to make sure that the fluid motion is restricted 
to a two-dimensional spherical surface. 

For such thin rotating shell, exact Rossby wave solu- 
tions^ of the fully nonlinear fluid equations are known from 
the geophysical literature [see, e.g., Silberman (1954)]. In 
Section II we review these solutions and study their proper- 
ties. We find that in our model Rossby waves carry zero 
physical angular momentum and positive energy. This is 
somewhat surprising, since Friedman and Morsink (1998) 
have shown that CANONICAL angular momentum and en- 
ergy of the r-modes are negative. It seems that one cannot 

* From here onward "Rossby waves" and "r-modes" will be used 
interchangeably. 
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equate canonical and physical quantities for r-modes; one 
can only do it if the Lagrangian perturbation of vorticity is 
zero [Friedman and Schutz, 1978(a)]. This is not the case for 
r-modes; this fact is intricately connected to the differential 
drift found by RezzoUa et al. (1999) and confirmed within 
our model. Our results do not argue against the presence of 
the CPS instability, since its derivation as given in Friedman 
and Schutz [1978(b)] relies only on canonical quantities. 

In Section III we switch on radiation reaction by assum- 
ing that the fluid is electrically charged. We derive coupled 
dynamical fluid — Maxwell field equations [see Eqs ( |l9| ) — 
(P5|)]. By studying the linear part of these equations, we 
then explicitly show the presence of the CFS instability for 
r-modes. We also show that the r-mode frequency is mod- 
ified by radiation reaction; the frequency shift is given by 
Auj ~ 7CFs(A/i?atar)'^'''^'^, whcrc I is the r-mode multipole 
order, 7cfs is the r-mode growth rate and A is the wave- 
length of the emitted gravitational wave [cf Eq. (^)]. 

We study the evolution of the star once the nonlinear 
coupling terms are included in the dynamical equations. We 
show that, to second order in the r-mode amplitude, radi- 
ation reaction slows down the star and causes it to rotate 
differentially, and that this is not related to the saturation of 
the r-mode. This is contrary to the phenomenological model 
introduced by Owen et al. (1998), and used extensively by 
Levin (1999). In that model it was assumed that, prior to 
saturation and in absence of viscosity, the loss of angular 
momentum to gravitational radiation is entirely accounted 
for by the r-mode growth. However, since (in our model) 
r-modes do not carry physical angular momentum, our find- 
ing only seems logical. The differential rotation is similar in 
magnitude, but different in origin from the differential drift 
found by RezzoUa et al. (1999). Both the slow-down and the 
differential rotation are described by Eq. ( p8[ ) of the text. 

Section IV provides a quick and somewhat superficial 
excursion into gravitational radiation reaction potential for 
mass current quadrupole moment, in the 3.5 post-Newtonian 
order. We use the formalism of Blanchet (1997) and others to 
show that in this order there is a full equivalence between the 
gravitational and electromagnetic radiation reaction. All of 
our results are thus expected to apply to the case of Rossby 
waves coupled to gravitational radiation, at least for I = 2. 
In particular, gravitational radiation reaction slows down 
the star and causes it to rotate differentially, prior to satu- 
ration of the r-mode growth. By carrying through results of 
Section III, we find that inertial frame dragging both from 
the background rotation and from the r-mode itself modifies 
the r-mode frequency by a factor of ~ iischwarzschiid/-Rstar, 
in qualitative agreement with results of Kojima (1998). 



2 ROSSBY WAVES IN A ROTATING SHELL: 
EXACT SOLUTIONS 

Thin rotating shells, and Rossby waves in them, have been 
studied extensively by geophysicists and meteorologists since 
the end of last century. In this section we follow closely the 
work by Silberman (1954); more original references can be 
found in that work. 

Consider a thin spherical sheet of incompressible fiuid 
of radius a, sandwiched between two spherical hard covers, 
so that motion of the fluid is restricted to a two-dimensional 



spherical surface. The fluid is rotating with angular fre- 
quency Q around the z-axis relative to an inertial observer. 
Since the fiuid motion is restricted to two dimensions, and 
the fluid is assumed to be incompressible, the fiuid velocity 
field in the co-rotating frame is completely determined by a 
stream function i/'(0, <j)) deflned on the fluid sphere: 



IdiP , 1 dip 

a do aamo do 



(1) 



where and Ve are the components of the fluid velocity 
along parallels and meridians, respectively. The vorticity of 
the fluid in the rotating frame of reference (called relative 
vorticity) is given by 



77 = (l/a)V^V, 

where is the Laplacian operator on a unit sphere: 



1 

sin ( 



d 



d' 



sin 6 d(j}'^ 



(2) 



(3) 



For inviscid, barotropiqlj flows, Euler's fiuid equations 
imply conservation of circulation: 



■ df= const, 



(4) 



C{t) 



where C{t) is a contour moving with the fiuid and iT'" is 
the velocity relative to an inertial observer. We derive the 
dynamical equations for an in incompressible spherical fiuid 
shell from this equation. In particular, Eq. ^ implies that 
the radial component of the absolute vorticity is conserved: 

^ (r? 4-20 cose*) = (5) 

[see, e.g., Longuet-Higgins (1964)]. Here d/dt = d/dt-\-v-V is 
the Lagrangian time derivative. We would like to stress that 
Eq. (jH]) is fully nonlinear, and the only assumptions made 
so far are those of incompressibility and zero viscosity. 

By expressing all velocity components and 77 in terms 
of the stream function [cf Eqs. (|^) and (^)], one can write 
Eq. (^) as a dynamical nonlinear equation for the stream 
function evolution: 

^vl^ + 2n^^ = ^{vli,,^}, (6) 

dt dcp am 

where ip = ip/a'^ is the reduced stream function, and the 
brackets are defined by {A,B} — deAd^B — deBd^A. The 
left-hand side of Eq. (^ contains only linear terms, while 
the right-hand side is the nonlinear advection term. 

Let us neglect this nonlinear part for a moment and 
look for a solution to the linearized equation in the form 



aQYime'''^* + cc. 



(7) 



where a is a constant. It is then easy to work out the dis- 
persion relation for the linear part of Eq. (|q) : 



2mSl 

i{i + iy 



(8) 



^ The exact condition for circulation conservation is Vpx Vp = 0, 
where p and p are the fluid density and pressure respectively. A 
fluid with a one-parameter equation of state, p = p(p), satisfles 
this condition automatically. 
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This is the well-known dispersion relation for Rossby waves. 
Moreover, since Yim is an eigenfunction of the Laplacian, 
the nonlinear term in Eq. (p|) equals zero when rp oc Yim- 
Therefore Eq. (Q) is an exact solution of the fluid equations 
of motion, for arbitrarily large a. 

Consider now the angular momentum and energy of the 
solution The total angular momentum of the fluid shell 
is given by 



L 



backgro' 



und H~ -^T 



(9) 



where -Lbackground is the angular momentum of the shell 
when the r-mode amplitude a is zero, and 



ii -modo = pa,^ / «</) sin^ 9d6d(j} 



(10) 



is the physical change in angular momentum of the shell due 
to the r-mode. Here p is the surface density of the fluid. It 
is then clear from Eqs (^) and (Q) that 

ir-modc ~ (11) 

for all m 7^ 0. Modes with m — have zero frequency and 
are called zonal currents; we will come back to them later 
when we discuss nonlinear evolution of the r-modes. On the 
other hand, the canonical angular momentum of a Rossby 
wave with I = m, computed from Eq. (3.4) of Owen et al. 
(1998), is 

^ [1(1 + l)f 4 2o no^ 

-t^canonical = Pd « (li) 

clearly different from the physical angular momentum. 
The total kinetic energy of the shell ^ is given by 



E = pa^ 



[(aSl sin 6 + v^)^ + d|] sin 



(13) 



For r-modes with m not equal to zero, the terms linear in 
V on the right-hand side of Eq. ( |l^ ) vanish after integration 
over (j>- Therefore, the total energy of the star with the r- 
mode is 



E -£^background -£^r — mo 

= pa^ I (aSl sin S)^ sin 



(14) 



+ pa 



The physical r-mode energy iJr-modc is therefore greater 
than zero. Both this and Lr-mode = [see Eq. ( |lll )] is con- 
trary to the common belief, which holds that iSr-modo and 
ir-modc are both negative. 

The source of misunderstanding is the confusion be- 
tween canonical and physical quantities, which are not equal 
to each other for r-modes. In their seminal work, Friedman 
and Schutz (1978a) have shown that canonical quantities are 
equal to the physical ones, so long as the Lagrangian change 
in vorticity is zero to second order in the perturbation am- 
plitude. The last condition cannot be true for r-modes; this 
fact is intricately connected to the work by RezzoUa, Lamb 



•t For an incompressible spherical fluid shell, both internal and 
potential energies are constant, and hence do not play any dy- 
namical role. 



and Shapiro (1999). These authors track the motion of a 
fluid particle for the case when the fluid stream function is 
given by t/) oc ¥22- They find that the particle experiences a 
drift along stellar latitude; the speed of the drift depends on 
the latitude. Therefore, the authors argue, the presence of 
an r-mode implies the presence of a differential drift in the 
star^. This is clearly incompatible with the zero Lagrangian 
vorticity perturbations; therefore, the Friedman-Schutz con- 
dition is not satisfied and one cannot equate canonical and 
physical energy and angular momentum|^ 



3 ELECTROMAGNETIC RADIATION 
REACTION FOR R-MODES 

Friedman and Schutz (1978b) had shown that the CFS insta- 
bility occurs whenever there is a mode with negative canon- 
ical angular momentum, which is dragged forward by the 
stellar rotation relative to an inertial observer, and which is 
coupled to some radiation field. This radiation field can be 
scalar, electromagnetic or gravitational — no matter what its 
nature, the CFS instability will be present. Papaloizou and 
Pringle (1978) investigated the CFS instability for f-modes 
coupled to a scalar field, and derived the growth rate of an 
unstable mode by introducing explicitly scalar-field radia- 
tion reaction. 

In this section we introduce electromagnetic radiation 
reaction (ERR) into our model. In the next section we will 
show that, for r-modes on a spherical shell, gravitational ra- 
diation reaction (in the lowest appropriate post-Newtonian 
order) and ERR are equivalent. Thus, all results derived in 
this section for the case of the ERR are applicable to the 
case of gravitational radiation reaction. 

We assume that the fluid is homogeneously electrically 
charged, with the charge q per unit massfj]. R-modes couple 
to electromagnetic radiation through time-varying current 



3 In private communications, RozzoUa, Lamb, and Shapiro had 
found some of relativity community to be sceptical about the re- 
ality of their claimed differential drift. This scepticism was due to 
the fact that RezzoUa, Lamb, and Shapiro derived the differential 
drift using the fluid velocities which were first-order quantities 
in the r-mode amplitude; however, the drift that they found was 
second-order in the r-mode amplitude. The sceptics thought that 
such procedure was inconsistent. However, for the case of a spher- 
ical shell, the fluid velocities of the exact Rossby wave solution 
are STRICTLY linear with respect to the r-mode amplitude, and 
are given by exactly the same expressions as used by RezzoUa, 
Lamb, and Shapiro. Therefore, at least for the spherical shell, the 
differential drift due to the r-mode is real. 

^ One other obvious example where there is no equivalence be- 
tween physical and canonical quantities is the shear sound wave 
in a solid elastic body. The physical motion of all the particles of 
the body is transverse to the direction of the wave propagation; 
therefore such wave has zero linear momentum along the direc- 
tion of propagation. Yet, the canonical linear momentum is not 
zero (in fact, it is p = Nhk, where A' is the number of phonons 
in the wave, and k is the wave vector). 

II Note that this situation is different from the usual MHD, where 
the current is due to relative motion of oppositely charged species, 
e.g., electrons and protons. Here, the charge current density is 
proportional to the mass current density. 
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multipole moments (since the fluid is assumed to be incom- 
pressible, the charge multipole moments are zero). We now 
derive dynamical equations of motion of the fluid shell cou- 
pled to an electromagnetic fleld. 

The circulation around a closed contour moving with 
the fluid is no longer conserved: 



dt 



^^-.dr=l '^.dr+ j . d ( ri. ] . (15) 

C(t) Jc(t) ™ Jc(t) ^ 

Here, as in Eq. (^), is the fluid velocity in the inertial 



-if 



frame of reference. The second term on the RHS of Eq. ( |15D 
is identically zero; the first one is in general not zero because 
of forces exerted on the fluid by the Maxwell field: 



d_ 

di 



■ df = 



C(t) 



q I E ■ df 

'c(t) 

q d 

r df , 

^ 'C{t) 



c 

B dA. 



X B-df 



C(t) 



(16) 



Here B ■ dA is the magnetic flux through the surface 

with boundary C, with dA being the area increment vec- 
tor. Equation (|l^) is the well-known Faraday's law, which 
states that the electromotive force (EMF) around a closed 
contour equals to the negative of the rate of change of the 
magnetic flux through the contour. We see that the presence 
of a Maxwell field modifies Eq. so that 



circulation + (g/c) ■ fiux = const. 
The analogue of Eq. ^ is then 

^ij) + 2Q.cose + ^Br) = 0, 



(17) 



(18) 



where Br is the radial component of the magnetic field in 
the inertial frame. Expressing velocity components in terms 
of the stream function derivatives [cf. Eq. (|l|)], we get 

at a<p c at sm6' 



(19) 



The angular operator is given by Eq. (^). Equation ( [l9| ) 
governs the fiuid motion under the action of an electromag- 
netic field, which, in turn, is generated by the fiuid currents. 

Let us consider a Rossby wave with a reduced stream 
function given by = + where 



(20) 



Then, as s hown in Appendix A [cf. Eqs (|Ag), (|A7|), (|A8D, 
(|AS|), and (AlC)], the radial component B,. of the magnetic 
field at the fiuid shell is 



Br — -Brot + Blm, 

where 

Brot = pi 2a cost 

3c 



(21) 



(22) 



is the radial component of a dipole magnetic field due to the 
uniform rotation of the charged fiuid shell, and 

Blm = xi'L + X*^Tm (23) 
is the radial component of the magnetic field produced by 



the Rossby wave itself. Here x = Xi + *X2, where, to lowest 
order in fea, we have 



Xi 



X2 



A-kI{1 + 1) qpa 

2r+i T' 



4nl{l + 1) qpa ,, .2i+i 
[(2Z + 1)!!]2 c ^ ' 



(24) 
(25) 



The wavevector k of the emitted electromagnetic waves is 
given by 

k ^ -{ioira - mQ.) = - — - - — . (26) 

c t t -f 1 c 

We now obtain the new dispersion relation for Rossby waves 
by substituting Eq. ( ^ ) into Eq. ( p^ ) and keeping terms 
which are linear in the mode amplitude a. The resulting 
linear equation is 



at ad) csmt 



Using Eqs (|2g), (g3|), (gj) and (|25), we derive the Rossby- 
wave dispersion relation [cf. Eq. (|8)]: 



2nm 



(1/3). 



e(fca)2'+V[(2; + l)!!]2' 



(28) 



(29) 



/(/ + 1) l + e/(2/ + l) 

where 

Anq^pa Q'^/a 

" C2 " Mc2 ■ 

Here Q and M are the charge and the mass of the fluid 
shell respectively; e is thus the ratio of the energy of elec- 
trostatic self-interaction and the rest-mass energy of the 
shell. For the case when a Rossby wave is coupled to grav- 
itational radiation, we will see that the analogue of e is 
~ -RschwarzschiidZ-Rstar iu the wcak gravity regime. For a neu- 
tron star, e ~ 0.4. Then, to flrst order in e, the angular 
frequency of the r-mode is 



where 



2nm 
1(1 + 1) 



1 



.7CFS 



2{l - 1) 
3{2l + 1) 



(30) 



(31) 



is the frequency shift of the r-mode due to the electromag- 
netic interaction, and 



7CFS ~ 



e(fca) 



21+1 



[(21 + 1) 



1112 



(32) 



is the growth rate of the r-mode due to the CFS instability. 
The angular frequency shift and the CFS growth rate are 
related by 



7CFS 



21 + 1 



(33) 



where A is the wavelength of the emitted radiation. 

Now we consider the terms which are of second order 
with respect to the r-mode amplitude a in the dynamical 
equation (^9|). The advection term (flrst term on RHS of 
Eq. [^) does not contribute to this order, but the second 
term on the RHS of the Eq. (|l|) does have a component 
which is proportional to a^: 



nonlinear term = 



csin^ 



(34) 
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{xi'L + X*4'lm, ^iln + i'lm} 



2ia O, Im{x)- 



■ {Ylm, i^im} ■ 



csin 9 

Let us focus, for concreteness, on the I = m — 2 mode, 
which has the largest growth rate (and hence is perhaps the 
most important astrophysically) . Then the nonlinear term 
in Eq. (|35|) becomes 



nonlinear term — 18a r27cFS 




(35) 



This term on the rig ht-hand side of Eq. (|l|) will act as a 
source for the left-hand side of this equation, thus creating 
a contribution to the stream function which is second order 



^(2) = Ki{t)Yw + K2(t)V30, 



(36) 



where, to leading order in e, ki and K2 satisfy the following 
evolution equations: 



dKi 
lit 

dK2 



9 



Q i27CFS, 



— a S27CFS- 

YTT 



(37) 



If the mode amplitude grows exponentially starting from a 
small value ao, i.e., a = aoe'''^^^', then by integrating Eq. 
( ^tI ) we get the following expression for the reduced stream 
function of the radiation reaction induced flow: 



(38) 



Let us discuss this equation. The first term on the right-hand 
side represents a uniform flow in the direction opposite to 
stellar rotation. It thus represents the spindown of the star. 
The angular momentum associated with this flow is 



(39) 



-^spindown 

in agreement with the canonical angular momentum of the 
r-mode itself [cf. Eq. (p^J. This is to be expected, as both 
^canonical and Lspindown must be cqual to the angular mo- 
mentum lost to radiation. 

The second term in Eq. ( ^ ) is a zonal current or, in 
another words, differential rotation. This term does not con- 
tribute to the angular momentum of the star. This differen- 
tial rotation is similar in magnitude to the differential drift 
reported by RezzoUa et al. (1999). However, its origin is 
completely different. In our case, the differential rotation 
is driven by the radiation reaction, whereas the differential 
drift of RezzoUa et al. is a kinematic property of the r-mode 
fluid motion, and is not at all related to radiation reaction. 

Note that both the uniform slow-down and the induced 
differential rotation are not related directly to the r-mode 
saturation, and are present well before saturation takes 
place. This is somewhat contrary to the phenomenological 
model of nonlinear behaviour of the r-mode instability by 
Owen et al. (1998), which was used extensively by Levin 
(1999). This model assumed that prior to r-mode satura- 
tion, the loss of angular momentum and energy carried off by 
gravitational waves was manifested by the r-mode growth, 
and that the background motion of the star was unchanged. 
The intuition for this model was based heavily on the belief 



that r-modes carry negative physical energy and angular mo- 
mentum. We now know that the latter is in general not true, 
so it should not be surprising that radiation reaction induces 
second-order changes in the background motion of the star, 
as well as drives the r-mode instability. Astrophysical impli- 
cations of this point are currently under investigation, and 
will be the subject of our next publication. 



4 GRAVITATIONAL RADIATION REACTION 
FOR L = 2 R-MODES 

An analogy between weak gravity and electromagnetism has 
been studied by many researchers [e.g., Braginsky, Caves 
and Thorne (1977), Thorne, Price and Macdonald (1986)]. 
Shapiro (1996) has shown that the Newtonian circulation 
around a closed contour comoving with a perfect fluid 
is not conserved in presence of a gravitomagnetic field; 
the conserved quantity (termed "relativistic circulation" by 
Shapiro) is a linear combination of the Newtonian circula- 
tion around the contour and a gravitomagnetic flux through 
the contour [cf. Eq. (4) of Shapiro (1996)]. This is very rem- 
iniscent of our conclusions for the circulation of a charged 
fluid in the presence of magnetic field; in fact, the two deriva- 
tions are almost identical. Since it is the circulation equation 
that determines the dynamics of our rotating shell, Shapiro's 
result is of great relevance to understanding this dynamics. 

In this section we use results of Blanchet (1997), Shapiro 
(1996), Asada et al. (1997), and RezzoUa et al. (1998) to in- 
vestigate the effect of gravitational radiation reaction on r- 
modes in a spherical fluid shell. We find that, for slow-motion 
systems, there is a close analogy with electromagnetic radi- 
ation reaction considered in the previous section. 

Following Asada et al. (1997) and RezzoUa et al. (1998), 
we consider a 3 -I- 1 splitting of spacetime. We write the 
square of the line element as 

ds^ = g^^dx'^dx'^ (40) 
= (a^ — Pi(3^)c'dt^ + 2Picdtdx^ + 'jijdx'dx-' , 

where a and are the lapse function and the shift vec- 
tor respectively, and yij are the spatial metric coefficients. 
For weakly gravitating (R ^ -Rschwarzschiid), slow-motion 
{v <C c) sources in which mass currents produce gravita- 
tional radiation, one can choose a gauge such that it is the 
time- varying shift vector /3 that plays a dynamically impor- 
tant role, relative to all other perturbations of the metric. 
For a periodic mass-current quadrupole moment, the shift 
vector consists of two parts, /3 = /^s™ -|- /^rr^ where 
1. The first part is the usual gravitomagnetic vector; its lead- 
ing term in u/c is given by 



c3 / r-r' 



(41) 



where v"" is the fiuid speed relative to an inertial observer, 
and a is the mass volume density of the source [cf. Eq. (3.4) 
of Blanchet (1997) when c — <x; Blanchet uses a vector 
potential V = -/3/4]. 

2. The second part is responsible for the radiation reaction; 
it changes sign under time reversal. Its leading term in v/c 
is given by [cf. Eq. (3.66) of Blanchet (1997) and Eq. (17.6) 
of RezzoUa et al. (1998)]: 
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16G 
45c8' 



:.(5) 



0' ^ /ii;„S ^»J'=-''*-^J'^fei ' 



where 



(fxeki(iXj)XkfJvr 



(42) 



(43) 



is the mass-current quadrupole moment. Here the super- 
script (n) stands for d" /df^, and the brackets around tenso- 
rial indices, (zj), indicate symmetrization over these indices 
(i.e. a(ij-) — l/2{aij + aji) for a tensor flij). 

The effective force F per unit mass that this shift metric 
perturbation exerts on the fluid is given by 



ot 



(44) 



cf. Eq. (12) of RezzoUa et al. (1998) and Eq. (1) of Shapiro 
(1996). Note that this expression for the gravitational force 
is equivalent the Lorentz force exerted by an electromagnetic 
field on moving charged fluid: 



q dA q T» „ 
-- — + -II'" X V X ^. 
c at c 



(45) 



In Equations ( |44[ ) and (^Bj) , the shift vector /3 is dynamically 
equivalent to c qA, where A is the electrodynamical vector 
potential. 

In our discussion of the motion of charged fluid on a 
spherical shell, we have shown that it was the radial compo- 
nent of the magnetic field that entered the dynamical equa- 
tions of the fluid motion. Likewise, an identical argument 
will work for the gravitational force given by Eq. (^). There- 
fore, the radial component of V x /3 enters the equations of 
motion of a gravitating fluid: 



06 at 



A^{V.V',^} + ^{b,Vi},(46) 



where b = (V x (3)r. 

Suppose that a single I = 2 Rossby wave is excited 
in a spherical shell with the surface mass density p, and 
that the wave's reduced stream function is given by ijjim = 
ar2y2me'"^™*. Then, as is shown in Appendix B, in the slow- 
motion approximation the radial component of the gravito- 
magnetic field generated by the fiuid motion is given by 



h = brot + b2m ; 
cf. Eq. (|l|). Here 

Orot = --Egrav— COSfcl, 

6 c 

and 

&2m = Xgrav'!/'2m, 

where, to lowest order in ka, 

T~t r \ 12ee:rav 

-Ke(Xgrav) = 

/m(Xgrav) = 

and 



(47) 

(48) 

(49) 

(50) 
(51) 

(52) 

Equations (^7|), (^sj), and (^) have the same structure as 
the analogous equations for the electromagnetic case, cf. Eqs 



3c 
225c 



{kaf, 



2GM Rs. 



chwarzschild 



(pi[), ([2^), and (pi^. By following the same steps as in the 
electromagnetic case, we work out the dispersion relation for 
I — 2 Rossby waves interacting with gravity: 



<^2m = 

where 



2n 

~3~ 



1 + 



.7CFS 

I 

<^2m 



(53) 



(54) 



is the shift of the r-mode frequency due to inertial frame 
dragging, which originates both from stellar rotation and 
from the mode itself; and 



7CFS 



c7m(Xgrav) 8eg,a 
-1^2m — 1^2n 



\ka\- 



6 225 ^^^^ 

is the growth rate of the r-mode due to the CPS instability 
driven by the gravitational radiation reaction. This growth 
rate agrees with the calculations of Lindblom et al. (1998) 
when one applies their Eq. (17) to the case of a massive 
spherical shell. 

The r-mode frequency shift due to inertial frame drag- 
ging was discovered by Kojima (1998). For the case of a real 
three-dimensional star, Kojima (1998) claims, and Beyer 
and Kokkotas (1999) confirm, that such shift causes the r- 
mode spectrum to be continuous. This claim, however, is 
not supported by calculations of Lockitch, Andersson and 
Friedman (1999). The issue of a continuous spectrum is not 
relevant for the spherical shell. 

Since the formalisms for a slow-motion gravitational ra- 
diation reaction from a mass-current quadrupole and for 
an electromagnetic radiation reaction from a charge-current 
quadrupole are identical in the structure of the dynamical 
equations, all of the conclusions from the previous section 
about the nonlinear electromagnetically-driven evolution of 
the I = 2 r-mode are also valid for a gravitationally driven 
/ — 2 r-mode, at least for the case when the r-modes are 
excited in a spherical shell. More specifically, Eq. ( |38| ) for 
the secondary reaction-induced flow is still valid. Therefore, 
gravitational radiation reaction will slow down the star and 
cause it to rotate differentially; the former will account for 
the loss of angular momentum to gravitational waves. Both 
the slow-down of the star and the reaction-induced differen- 
tial rotation are not related to the nonlinear saturation of 
the r-mode growth. 



5 CONCLUSIONS AND BETS 

This paper has studied the issues of principle for the r-mode 
instability in the nonlinear regime. Although we considered 
only a special case of r-modes excited in a spherical ro- 
tating shell, we bet that most of the lessons learned from 
the simple model will apply to real stars. In particular, r- 
modes in general do not carry negative physical angular mo- 
mentum and energy; the radiation reaction causes the star 
to slow down and rotate differentially prior to the r-mode 
saturation. Therefore, phenomenological nonlinear evolution 
Equations (3.14) — (3.17) of Owen et al. (1998) need to be 
reconsidered. 

We believe thus that our conclusions argue in favor of 
Spruit's conjecture (1999) that the r-mode instability causes 
differential rotation in the star. Spruit has modeled differen- 
tial rotation as relative motion of two spherical shells, while 
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we find difTerential rotation in the lateral direction. We be- 
lieve that, in the three-dimensional case, both radial and 
lateral differential rotation will develop. A more detailed 
discussion of this and other astrophysical consequences of 
our formalism is a subject of a future publication. 

We want to thank Luciano RezzoUa, Frederick Lamb, 
and Stuart Shapiro for showing us the advance draft of their 
manuscript, and for detailed comments on our first draft. 
Discussions with John Friedman, Peter Goldreich, Keith 
Lockitch, and Kip Thorne have been very helpful. Alan 
Wiseman has pointed us to useful references on radiation re- 
action. YL was supported by the NSF grant AST-9731698, 
and by the Theoretical Astrophysics Center at UCB; GU is 
a Fannie and John Hertz Foundation Fellow. 



REFERENCES 

Andersson, N. 1998, ApJ, 502, 708. 

Andersson, N., Kokkotas, K. D. & Schutz, B. F. 1999, ApJ, 516, 
307. 

Asada, H., Shibata, M., & Futamase, T., Prog. Theor. Phys. 

Suppl., 96, 81. 
Beyer, H. R., Kokkotas, K. D. 1999, MNRAS 308, 745. 
Blanchet, L. 1997, Phys. Rev. D, 55, 714. 

Braginsky, V. B., Caves, C. M., & Thorne, K. S. 1977, Phys. Rev. 
D, 2047. 

Friedman, J. F., & Morsink, S. M. 1998, ApJ, 502, 714. 
Friedman, J. F., & Schutz, B. F. 1978a, ApJ, 221, 937. 
Friedman, J. F., & Schutz, B. F. 1978b, ApJ, 222, 281. 
Kojima, Y. 1998, MNRAS, 293, 423. 
Levin, Y. 1999, ApJ, 517, 328. 

Lindblom, L., Owen, B. J., & Morsink, S. M. 1998, Phys. Rev. 
Lett., 80, 4843. 

Lockitch, K. H., Andersson, N., & Friedman, J. F. 1999, unpub- 
lished work; preli minary results are reported in K. H. Lock- 
itch's PhD thesis, |gr-qc/990902£ . 



Longuet-Higgins, M. S. 1964, Proc. Roy. Soc. London A, 279, 446. 
Owen, B. J., Lindblom, L., Cutler, C, Schutz, B. F., Vecchio, A., 

& Andersson, N. 1998, Phys. Rev. D, 58, 084020 
Papaloizou, J., & Pringle, J. E. 1978, MNRAS, 184, 501. 
Rezzolla. L.. Lamb. F. K., & Shapiro, S. L. 1999, submitted to 

ApJ, istro-ph/9911188 
Rezzolla, 



L., Shibata, M., Asa da. H.. Baumg arte. T.W., 
Shapiro, S. to appear in ApJ (gr-qc/9905057). 



Shapiro, S. L. 1996, Phys. Rev. Lett., 77, 4487. 
Silberman, L 1954, Journal of Meteorology, 11, 27. 
Spruit, H. C. 1999, A&A, 341, LI. 

Thorne, K. S., Price, R. H., & Macdonald, D. A. 1986, Black 
Holes: the Membrane Paradigm (Yale University Press, New 
Haven), Chapter 4. 



APPENDIX A: MAGNETIC FIELD 
GENERATED BY A ROSSBY WAVE IN A 
CHARGED FLUID. 

In this Appendix we find the radial component of the mag- 
netic field generated by a Rossby wave in a charged fiuid. 
We thus derive Eqs (|25|) of the text. 

Suppose that a single Rossby wave is excited in a ro- 
tating shell, and that its reduced stream function is given 
by 



The radial component of the magnetic field produced by 
such fluid motion can be found by using the multipole for- 
malism discussed in Sec. 16.5 of Jackson (1975). In particu- 
lar, using Eq. (16.87) of this reference, we find: 



Br = Brot + Bi 

where 



-Orot — 



and 



R -^P 



20, cos 6 2 . 
—a am 



~ik\ r — r' \ 



V X va^ sin O'dO'o 



(A2) 



(A3) 



(A4) 



Here Brot is the dipole magnetic field due to the uniform 
rotation of the charged shell, while Bim is the field due to 
the Rossby wave. The wavenumber k is that of the emitted 
electromagnetic radiation; 



mO,). 



(A5) 



Equations (A3) and (A4) are evaluated by noting that 



V X i7= aVa-lplni 
— ik\ r — r' I 



al{l + + C.C., and that 

Yim{e', 4,') sin e'de'd(t)' (A6) 

r — r' I 

^ 47rikh'i'\ka)ji{ka)Yi^{e,cl)), 



where ji and hj are spherical Bessel and Hankel func- 
tions respectively (see e.g. Eqs (16.9) and (16.10) of Jackson 
(1975)). After some algebra, we get 



Brot ~ ^^pHacosS, 
3c 

and 

Blm = XV'iln +X>im 

Here x — Xi + *X2, where, to lowest order in ka, 
47vl{l + 1) qpa 



Xi = 
X2 = 



21 + 1 c 
4tvI{1 + 1) qpa 

[{2rTTw~ 



(ka 



.21+1 



(A7) 
(A8) 

(A9) 
(AlO) 



(Al) 



Thus, we have derived Equations (^), (p^), (^^, (p^, and 
(H) of the text. 



APPENDIX B: GRAVITOMAGNETIC FIELD 
GENERATED BY AN i = 2 ROSSBY WAVE IN 
A THIN SHELL 

In this Appendix we derive Eqs (^, (^), (|^, and ^ 
for the radial component of the gravitomagnetic field, b = 
(V X /3)r, generated by a rotating massive shell in which an 
/ = 2 Rossby wave is excited. 

Evaluation of the part of b which is not responsible for 
radiation reaction, to leading order in v/c, is straightfor- 
ward. The relevant part of the shift vector, /3nonradiativo, is 
given by Eq. (^l|) of the text, which, up to a constant fac- 
tor, is same as the nonradiative part of the electromagnetic 
vector potential; 
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A 



1 



nonradiativc 



Therefore, 

/5non radiative — 
^nonradiativc 



AG J 

o -^nonradiativc ; 

c^q 

J^rnonradiativc ■ 

c^q 



(Bl) 



In the new coordinate system we can now evaluate b{r): 



(B12) 



(B2) 




cos U\ 



r-V ^ al,„Y2m{9i,(t)i) 

m=-2 



(5) 



sin 9id9ic 



Here, as in the text, u is the mass volume density of the 
fluid (trivial changes must be made for the case of a two- 
dimensional sphere), and q is the charge per unit mass of the 
fluid. By using Eqs. @, (||), (H), and (|^) of the text, 
and substituting I — 2, we get 



In Eq. (BI3), all terms with non-zero m vanish after inte- 
gration over (p; only the term with m — contributes to the 
integral. Integrating the remaining term, we get 



bir) 



48Gp 8 4/51 



(5) 



(B13) 



___4 n 12£grav 

t'nonradiativc — ^ f^grav COS u -f- 

3 c 3c 



2GM/{c^a) = Rsc 



V'2m, 
tiild/ffl- 



(B3) 



Substituting Eq. (Bll) into Eq. (B13), and evaluating the 
fifth time derivative, we get 



where egr 

Now we shall derive the expression for the part of the 
radial gravitomagnetic field which is responsible for radia- 
tion reaction; we shall denote it by The radiative part of 
the shift vector is given by Eq. (U2) of the text, which can 
be rewritten as follows; 



(B14) 



Here egrav = 2GM/(c^a) = -Rschwarzschiid/a. Equations (|B3[) 
and (|B14|) are equivalent to Equations (^), (p8|), (p9|), (^OD, 
and ([51[) of the text. 



or, in index- free form 



(5) 



(B4) 



where 

j' = r X v{r) = —'V'ip{r). 



(B6) 



In this expression, vectors r, r*, and the fifth time derivative 
are defined relative to an inertial observer at rest. The radial 
component of the radiation reaction gravitomagnetic field is 
given by 



48G'cr 



6" = -(r-Vx/?^ 

a 45c*a 

or, for the case of a spherical shell, 
48Gpa 



a{f- r)[f- Vip(f')]d^x' 



(5) 



.(B7) 



45c8 



(r-r)[f-Vi/)(r)] sin e'dO'i 



(5) 



(B8) 



Let the stream function tp he sl linear combination of 
I = 2 spherical harmonics: 



2 



lP{9,<l>) = a2my2m(6l,(^). 



(B9) 



Let us choose a coordinate system Xi,yi,zi (with polar an- 
gles 9i and (pi), such that the axis zi is along r. In this 
system, is also a linear of Z = 2 spherical harmonics, but 
with different weight coefficients a2m'- 



■^{n) = ^ al^Y2m{9i,(pi 
For r 1 = r (i.e. for 9i — 0), 

5 

ip{r) = 020^20(0, (pi) 



(BIO) 



(Bll) 
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